{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "name": "ThermalStimulation.ipynb",
      "provenance": [],
      "authorship_tag": "ABX9TyNa7LYgWsaPY3bqss7sZkuP",
      "include_colab_link": true
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/Divyanshu-ISM/Oil-and-Gas-data-analysis/blob/master/ThermalStimulation.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "7pukrCasp_T4",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "import numpy as np\n",
        "import matplotlib.pyplot as plt\n",
        "\n",
        "%matplotlib inline"
      ],
      "execution_count": 45,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "Vo-nujNUqwMy",
        "colab_type": "text"
      },
      "source": [
        "# Radial form of Darcy's Law\n",
        "\n",
        "#1. $ q = \\frac {kh (Pe - Pwf)}{141.2 \\mu B ln(\\frac{re}{rw})} $\n",
        "\n",
        "#2. $ \\triangle P  = \\frac {141.2q \\mu B ln(re/rw)} {k h} $\n",
        "\n",
        "#3. $ Pr = Pwf + \\frac {141.2q \\mu B ln(r/rw)} {k h} $\n",
        "\n",
        "#4. $ q = \\frac {kh (Pe - P)}{141.2 \\mu B ln(\\frac{re}{r})} $\n",
        "\n",
        "#5. $ \\triangle P  = \\frac {141.2q \\mu B ln(re/r)} {k h} $\n",
        "\n",
        "#6. $ Pr = Pe -  \\frac {141.2q \\mu B ln(re/r)} {k h} $\n",
        "\n",
        "\n",
        "\n",
        "\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "g4ZrxkT-qsRr",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "###### CODE As follows #######"
      ],
      "execution_count": 64,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "fBiBvXe7vIGk",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "# See that the blue curve is of the unstimulated-thick oil scene. You can see that oil's pressure diminishes even before \n",
        "# it reaches the well. Where as, the orange curve (hot oil) is still energetic by the time it reaches the well "
      ],
      "execution_count": 65,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "xd9d-uXIvWEA",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 421
        },
        "outputId": "334f6cb3-ba16-4c95-fee8-9c1f6bb070e8"
      },
      "source": [
        "plt.figure(figsize=(10,6))\n",
        "\n",
        "plt.plot(r,Pr_normal,linewidth='3',label=' without stimulation, oil viscocity = 1000 cp')\n",
        "\n",
        "plt.plot(r,Pr_stim,linewidth='3',label='with stimulation, oil viscocity = 10 cp')\n",
        "\n",
        "\n",
        "plt.grid()\n",
        "\n",
        "\n",
        "plt.legend(loc='best')\n",
        "\n",
        "plt.ylim(0,3000)\n",
        "\n",
        "plt.xlim(0,1000)\n",
        "\n",
        "plt.ylabel('Pressure (Psia)')\n",
        "\n",
        "plt.xlabel('Radius from the wellbore,r(feet)')\n",
        "\n",
        "plt.title('Qualitative Effect of Thermal Stimulation : Magic of oil viscosity reduction!')\n"
      ],
      "execution_count": 63,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "Text(0.5, 1.0, 'Qualitative Effect of Thermal Stimulation : Magic of oil viscosity reduction!')"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 63
        },
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAncAAAGDCAYAAABJITbwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdd3zddb348dc7u9k7TZq26UwXHbTQMi2WKcvBXg4UvMoFFfUq199VRLx4VcSriCCgrILIuAxRhhAKtLR0r3Q3aWaz904+vz8+3yQnyUlykubknCTv5+NxHuec7zrvc75nvM9nijEGpZRSSik1PgT4OgCllFJKKTVyNLlTSimllBpHNLlTSimllBpHNLlTSimllBpHNLlTSimllBpHNLlTSimllBpHNLlTfYjIl0TkQ5f7dSIycwSP/0cR+X8jdbwhPO7PRKRMRIqd+58TkTzn+S0b7Xhc4hp2HCLyFxH5mbdiOxEjGZuITHNen8CROF6vY/d4vw9j/3+IyBdHMiZ/ISLXi8hbXjhuj8/iEPft+j4a7ntMRM4Skf1D3e9Eeev1HC4RyRGRc71wXJ98x6tumtyNUc4P0i4RaRCRYhH5g4jEeOOxjDGRxpgjzuMO6cvU3Q+nMebrxph7RjpO54uq0fny77z83lk3DbgTWGCMmezs8ivgNuf5bTuBxzUiMvsEQncbh0tC03kxIlLvcv+sE3hMnxKREBH5tYjkO88lR0QecFnf40fHGHPMeX3afRNxV1w/EZGnXZcZYy4yxjzhg1iynPfEkl7LX3aWrz7RxzDGPGOMOf9Ej+Oqn8/iUGLq+j4aLmPMB8aYzBM5xjAft8frOQLfHT7n7e94533+1ZE41kSiyd0YJCJ3Ar8AvgfEAKuADOAtEQn2YWj+4FLny7/zcpuzfBpQbowpcdl2OrBn9EPsw20cLglNpDEm0lm8xGXZByMZhDdKxQbwQ2AFcCoQBawGto7i448XB4CbOu+ISAJwGlDqs4gG5+6zqE6QiAT5OgblPzS5G2NEJBq4G/h3Y8w/jTGtxpgc4CpgJnCds12PEjYRWS0i+S73fyAih0WkVkT2isjnBnhMIyKzReQW4Hrg+05py2sDHUtE5gN/BE5ztq/qHZuIZIvIJS6PFSQipSJysnN/lYisF5EqEdkxnNIIpwTobSDNieNZEakDAoEdInLY2S5NRF50Hv+oiNzucoxAEbnL5XluEZGpIrLO2WSHc+yr3Tx+gIj8SERyRaRERJ4UkRgRCXUXxzDEicjfnbg2isgsl8eeJyJvi0iFiOwXkatc1v1FRB4SkTdEpB44xykx+56I7HRKCR8TkRSxVY+1IvKOiMS5HONvYkuOq0VknYgs9DDmU4CXjTGFxsoxxjzpHPMpbALwmvOafl9EMpz3YZCzTZbYqr31ne9FEUkQkWdEpEZEPhGRDGfbHvu67O+2NEBEfiu2mrzGOc9nOcsvBO4CrnYec0fvY/V3rnvF8UUROSa2WvI/PXy9+vOME09nYn4t8DLQ4vJ8ThWRDc5nqEhEfi8iIS7rz3feG9ViawDed3k+vZtoLHR5Px0Xkbv6eQ1jnOde6rwWP3Jem96fxb/0s//XROSQ8zivikiay7pBS7ucz1aViCxyWZYktmQ/Wfp+H/6HiBQ47/H9IrLGWe72c++sO915n1U716e7HO9LInLE2eeoiFzf+/UUN98dIrJbRC51OU6w8z7p01yj8zk4sRcDf3Ze487v43IReV5E4l32udE5H+W933sy+G/GVBF5yTmn5c77aNDveA/P59dF5KBzzh4UERno/CoPGGP0MoYuwIVAGxDkZt0TwDPO7b8AP3NZtxrId7l/JZCGTfCvBuqBVGfdl4APXbY1wGx3xx3qsXofA/ivzpid+xcD2c7tKUA58Bnn2Oc595P6eW1ygHP7Wdfj+bt5XgHAFieeEGyifAS4wFn/PWAXkAkIsARI6H2cfh77K8Ah55iRwEvAU+7iGOTc99nOeS3LsSVgQdgf++ecdRFAHvBlZ90yoAxbHda5bzVwhvP8w5zX8GMgxXn9S7Alasuc9e8CP+713KKAUOABYLu78+zmufwIOAZ8AzgJkIHOJbZk2uC874Es5zWdhS293ostxTrXea5PAn92t6/L/l/t5/1+A5DgHOdOoBgIc9b9BHi6V6yux+r3XLvE8SdgEvY91AzM7+c1ug7YOcD7IQv4KvAWcJGzbBO25C4fWO0sW44t3Q9yYsgGvuWsSwRqgM876+8AWt29Ns55LnJekzDn/sp+YnsSeMXZJsM5Nzf391nste+nse/Tk7Hvq98B6zz9PnLZ7nHgXpf73wT+2TsG7Gc6D0hzOU+zBvrcA/FAJXCj87pd69xPwH7uaoBM5xipwMLBvlud+98H/upy/3Jg1wDfaW3YWpxQ5z11B/bzm+4sexh41tl+AVAHnO2su9/Z/1x3r2Wv1ygQ2AH8xnl+YcCZ7p5T72N5eD5fB2Kxf+pKgQvdfb704vlFS+7GnkSgzBjT5mZdEZDkyUGMMX8zttSkwxjzV+AgNkEYshM81lrgMhEJd+5fBzzr3L4BeMMY84Zz7LeBzdhkrz//5/z767x8zcM4TsEmjT81xrQY26bnT8A1zvqvAj8yxuw31g5jTLmHx74euN8Yc8QYU4etkrxGRq4a5WVjzCbnPfEMsNRZfgmQY4z5szGmzdj2fC9ik/FOrxhjPnJe3yZn2e+MMceNMQXAB8BGY8w2Z/3L2EQPAGPM48aYWmNMMzbxWSKetf38b+yP0vXYc1ogQ++U8GdjzGFjTDXwD+CwMeYd53X4m2ucQ2GMedoYU+68Zr/G/iB52j7Lk3N9tzGm0RizA/uDucTdgYwxa40xiz14zCeBm0RkHhBrjNnQ6zhbjDEfO88nB/uD/yln9WeAPcaYl5zX7X+xyaw7lwDFxphfG2OanPO+sfdGTiniNcAPnW1ygF9jEyFPXA88bozZ6ryvfogtGcrwcP9Oa+n+/IL9blnrZrt27DleICLBxpYid5ai9/e5vxg4aIx5ynldnwX2AZ2lbh3AIhGZZIwpMsZ42vzjaeAzYmtowL5mTw2wfQf2z1azMaYR+Drwn8aYfJfP5BXO++8K4HVjzDpn3f9z9vfEqdg/8N8zxtQ759/TTkienM/7jDFVxphjwHt0f4epYdLkbuwpAxL7SQxSnfWDEpGbRGR7ZxIELMImjkN2IscyxhzCliRc6iR4l9H9BTwduNI1WQPOxD7P/nzWGBPrcvmTh09jOraqyPWx7sKWYAFMBYZbbZoG5Lrcz8X+209xv/mQuf4YN2BLjMA+p5W9ntP1gGsj9jw3xzvucrvRzf1I6Kqyus+pAqrBlraBB+feGNNujHnQGHMG9h/7vcDjTjWPpzyKc6hE5LtimwtUO69ZDJ5/Njw51/2dr+F6CVs6chtuEgERmSsir4utPq8Bfk7380nD5T1gjDHYUj93PP0MJALB9H0dpniwb2dMXfs6SXL5EPbv9B4QLiIrnURiKfbPSQ/Od9C3sIlQiYg851Jt2N9z7n2ece5PMcbUY2swvg4UiW0yMc+TgI0xhcBHwBdEJBa4CPuHrT+lLn/KwH7mX3b5vGdjk9cU+p7reuzr6ompQG4/hQqD8eR8jvRnYsLT5G7s2YCtyvm860IRicR+EWQ5i+qBcJdNJrtsOx1bKnUbtmoxFtiNrXYYjOn1uIMdyzC4Z7HVGpcDe50vW7BfRE/1StYijDH3eXDMocoDjvZ6rChjzGdc1s8aYP+BFGK/dDtNw1aHHHe/+YjJA97v9ZwijTH/5rKNJ+enP9dhz9m52AQow1k+pPYyTinWg9hqrQUjEFdv9c6128+DK7Ht676PbcMa57yfq/H8/Tzq59oY04Atufw33JfyPIQtVZpjjInG/mnpfD5F2Co8AJy2Tul9jmDlYaubB1OGrdrt/ToUeLAv9HoNRSQCW93p6f6A/QMBPI/9brkWW2pV28+2a40xZzqPa7ClytD/5773eQaX52iMedMYcx72j+g+7Hekp57A1lpcCWxwStD70/v9mIetonf9zIc5xyjCJmkAOH+mE1z27fc3wznutH4KFYb0mRjq+TTGrDbGPOrJtqqbJndjjFMFdTfwOxG50Glwm4H9Eiuj+1/edmzxfryITMb+M+0Ugf1AlgKIyJexpW2eOE7PL/jBjnUcSBeXBtxuPAecj/1xcq02eRpboneBU0oU5jTy7e/H50RsAmqdxsmTnMdbJCKnOOsfBe4RkTliLRbbMxH6via9PQt8W0RmOEn4z7HtaobzL3goXgfmOo2og53LKUMsHRtIFPaPRjn2R+Hnnu4oIt9yzuUksZ1ovugcr3MomMFeU48ZY0qxPyQ3OOf1K/SfqEdhk7FSIEhE/guIdll/HMgQkf6+O311ru8CPuVUgfYWhW0DVueUILkm938HThKRzzo/3N+kn8QX+35Kdc5dqIhEicjK3hu5JFX3OttMB76D/Tx74lngyyKyVERCsa/hxn6e22DWYkvRrsd9lSwikikin3Yeqwlb6ttZXdnf5/4N7GfrOuf9ezX2j8nrYjsgXe4kMc3Ydm79VX+6e5//H7Z92h3YKveh+CP2dZ/uPLckEbncWfcCcImInOl8H/+UnjnAQL8Zm7DJ4X0iEuF8F5/h8hwG+o4fyfOpPKTJ3RhkjPkf7Jf5r4Ba4Cj2x/Vcp6gd7D/4HdiqsreAv7rsvxfbBmYD9oN5ErYqwBOPYdumVInI/3lwrHexw3wUi4jbKmNjTJGz/+m94szDlgzdhf2xzcM2cB7ofdvZw7Lz0qcapp8Y2rFtipZiX88y7Bd7Z/ux+7E/WG9hfygfwzZgBlud84TzmlxFX49jz8c659hNwL97EteJcEopzse2OyrEVn10Nr4eCU9iq1sKsB0aPh7Cvg3Y900x9rX+JvAF0z1+2X8DP3Je0++OQKxfw753yoGFwPp+tnsT+Ce2A0Au9ly5Vl3/zbkuFxF3Q7eM2LkWO+CtR221jG3z2l8bqO9iS1lrsSVIrp+xMmwJ0f9gX5sF2DaQzW4eoxbbqelS7Hk7CJzTz2P+O7Yk6AjwITaxetzD5/IOtj3Yi9iEYhY92855zGkTWI+tGvxHP5uFAvdh34fFQDK2XRj087l32t1dgu1cUo4t7b3EeT0DsMlsIVCBbd/omlC7+gm9vjuctnMvAjOwVe5D8VvgVeywWLXYz+RK57h7sJ+ztdjXtZKeVfAD/Wa0Y8/7bGxHqHxs0gyDfMef6PmUcTxIuDeJbWKhxjKntOynwBlOg1SllBoyp0QyH7jeGPOer+OZqJwS47nGmBt8HYsam3TQw3HAGPNnEWnDlnxpcqeU8piIXABsxFZHfg/bHm8opbBqBIkdl+5mPO9drFQfXquWderkN4kdeHaPiNztLJ8hdqDVQyLy1856eqcNx1+d5RvFpZu0iPzQWb7f+SJSvThd8p/zdRxKqTHnNGyP0DJs1dtnnapBNcrEDt2UB/zDGLNusO2V6o/XqmVFRIAIY0yd2CmxPsQ2EP0O8JIx5jkR+SOwwxjzkIh8A1hsjPm6iFwDfM4Yc7WILMA2yOwcZ+cdbHG1T+eYVEoppZTyR14ruXMGfKxz7gY7F4Mdj+kFZ/kTwGed25c793HWr3ESxMuxI+43G2OOYkd/H9Zgu0oppZRS451Xe8s6ww5sx05h9Da26L/KZViAfLoHMpyC0yvNWV+NHQuna7mbfZRSSimllAuvdqhwqk6Xih1p+2XAo1G6h0PspPa3AISFhS2fNm2atx5KeVlHRwcBATpKz1il52/s0nM3tun5G5yYdsR0ONftg9y31yM7prp7W4o6yowxHk0f6olR6S1rjKkSkfewDXdjRSTIKZ1Lp3uU6gLs6Nn5zmCaMdjxgzqXd3Ldx/UxHgEeAcjMzDT79+/31tNRXpaVlcXq1at9HYYaJj1/Y5eeu7FtQp2/jnZorIKGcmissNddlwrn0mtdYxWDJ2oCBDqXEwjPCDWEU2kiqSaSKhNJFRFUmiiqicCExRERm0RcQjIpKal8as3FvaezOyFeS+5EJAlodRK7SdjBL3+Bne/vCuysBF8EXnF2edW5v8FZ/64xxojIq8BaEbkf26FiDna0bKWUUkqNdcZASx3Ul3UnYvVl0ODcr3eTwHmUqI2MtsBJ1AdGU2kiKWkLp6QtnCoTSSVR9tpEUuUkcJVEUWkiqSGCDgKICg0ic3IUmZOjmJcazaLJUcxNiSJmUrBXY/ZmyV0qduTtQGzbvueNMa+LyF7gORH5GXaqocec7R8DnhKRQ9hRva8BO6q2iDyPHQG/Dfim9pRVSiml/FRHe3fJWUOZm6StvGfi1lAO7X0mRfGOSXEwKR7CE5xLPEyKoyk4loLmSeQ0hHKwNojdlYHsKA+kpC2cZgaaPbNbRkI481OjXS5RTImdhO0bOrq8ltwZY3YCy9wsP4Kb3q7GmCbsNDjujnUvcO9Ix6iUUkqpQXQma/WlTrJWapOy+lKXZS5JXGMlo1KqFhbbM0lzve6dwIUnQFgsJiCQgqpG9hTWsKegmr1FNWTvq6WgyvOhHScFBzIvNaoriVuQGkXm5GgiQ/1nXgj/iUQppZRS3tfRAU1VTkLmJkGrL+2uFq0vtYmdt5O14HAIT7SJWESik5QlQkRC9+1wl4QtLBYCB05h2jsMR8vq2FNQw57CUvYUHmJPYQ1VDa0eh5UWE9anNG56QgSBAaNfGjcUmtwppZRSY11bCyHN5VC43SZkdSXdiVtdCdSXQF1nIlcO3m7dNCnOScgSXJK1ztudyxO6b4eEn9DDNbe1c/B4HXsKq9ldUMOewmqyi2ppbPXseYYEBjAnJbJHErcgNZrYcM+qZP2NJndKKaWUP2ppsElZfVmvBM1J3Dpv15VAUxWng+2S6A2T4iAiySlN67wkOcsSnNtO4jYpbtBStRNR19xGdpGtVt1TWMPuwhoOldTS2u5Z6WJUWBAL06JZmBbDwrRoFqRFMyspkuDA8TOMjCZ3Siml1Ghpa+lOyOpKoO64y7XL7fpS24PUW0Kju5OxzsTMbcKWZKtDA73bu7M/VQ0t7OpM4gqq2VtYw9HyejydOTUpKpRFLoncoikxpMf5ppPDaNLkTimllDoRHR12qI7eCZq7pK2x0jsxSAAtQVGExKVDZBJEJNvErM9t5xIU6p04TkB1Qyu7C6vZmV/N7oJqdhZUkVfheUeHafHhXQncgrRoFqZFkxwV5sWI/Zcmd0oppZQ77a02Iasttpe6YpfbrklbiXfasAWGOIlZIkQm973dlawlQ3g869d9MGYGMa5pamV3QTW78qvZVWAvueUNHu0bIDA7OZJFaZ1JnL329thxY4kmd0oppSaWtubuBK22CGqP903caotsx4ORJgFOYpYMkSnOJbnvdUQShMXAOKg+rG1qZU9hTY9E7mhZvUf7hgQGkDk5ikVTuqtV502OIiz4xGaQGO80uVNKKTU+tLc6CVoR1BR23+6dxHmjajQstleC5i5pS7Ht1wLGb2JS39xmE7mCanblV7HTSeQ8aSMXHChkTo7ipCmxnDQlhsXpMcxNiSIkaPx0dBgtmtwppZTyb8bYcdlqiqC20LkudrntXNeXMqLjsUmALUGLmgyRkyEqBaJSbZLWuSzSKYXzwzZs3tbc1s6+olp25FexPa+KXfnVHCqt8yiRCwoQ5qZEsTg9hkVOIpc5OYrQoPGb+I4mTe6UUkr5TnurU6pW7JS2FblcuyRubZ43rB+UBDoJWq9krXcSF57o1SE9xpKODsPR8np25FWxI6+K7fnVZBfW0NLeMei+gQHCnORIFqfHcNKUGE5Kj9WqVS/Td61SSinvaG9zkrQCe6kusIlbTX737brjjFxpm9hStKhUiE5zErY0m6xFTu5O4MITxnXV6EgoqWlie14VO/Kr2JFXzY78Kmqb2gbdL0BgTnIUJ3UlcjEsSI3WRG6UaXKnlFJq6NrbbGJWUwDV+U7S5prEFdj1ZvCSHY8ER0B0at/ELTq1+zoyxWfjsY1ldc1t7OxM4pyErqi6yaN9MxLCWTI1liXpsSxOt71Ww0M0tfA1PQNKKaV6MsbOJ1qd51zyuy+dSVxt8QgN/yE2KXNN0twlbqHR46LnqK+1tnewv7jWlso5idzBEs/aySVEhHQlckumxrAkPZa4iLE5Pdd4p8mdUkpNNO1tti1bdT5U5UH1MajOZ/GRHbC7wSZ0rZ6NOTaoyBRb0hY9BWLS+96OStXSNi8qrm5i67FKtuZWsi2vit0F1TS3DV6aOik4kJOmxNgkzknoJsLMDuOFJndKKTXeNNc5JW1OyVtVXvf9qjyb2LmpLo0f6uNEJNlELXoKxEzpdTvNlroFacnOaGnrMGw7VsnWY1VsPVbJttxKCj2oXg0QmJsSxdKpsSyZGsvSqbHMSY4kaBzNtTrRaHKnlFJjTUs9VB2Dylx7XZULlTn2dnXeyIzjFhIJMVMhdqotZYtJh+j07iQuKhWCJ+bUTv7CtVRu67FKduY30PbW+kH3S4+bZJO4dJvMLZqi7eTGGz2bSinlb9panFK2XCeBy3VJ5nKd8dxOUGSKk7R1JnBT2ZVXxUlnXGTvh8VqGzc/0tLWwZ7C6iGXyk0KDmRxegwnT4/j5GlxLJ0aS1LUxBuTb6LR5E4ppUZbR7sdIqQzWetxfazfalOPBYbY0jUnabOX9O770VPclrqVN2ZB6uLhP64aMb1L5XYX1tDiQVu56QnhLJsa25XMzZscpdWrE5Amd0op5Q2tTbaqtPIoVByBCue68qht99bROvxjBwQ5ydp0iJvuXGdA7DR7iUiGAP1BHyvaOwz7i2vZnFvBJzmVbMmpGFapXGPeXi674JxRiFj5O03ulFJquJpq3CRvOfa6ppDhD84rtkNC7LReCZxzHZ2mg/COYY0t7WzPq2KLk8xtza2ktnnwAYKnJ4Rz8rQ4Tp4WyzI3pXJZJdneDFuNIZrcKaVUf4yBhvK+JW+d9xvKhn/s8EQnWeudwGXYUrkJOFfpeFVe18zm3Eo259hkbndBNW0dAyf+vUvllk2LJTFS3xPKM5rcKaVUYxWUH4byg1B+CMoOdidwLbXDO6YE2PZt8TMgfibEzbC342bYBC40ckSfgvIPxhhyyxv4JKeCzTmVfJJbwZHS+kH3S44K5ZSMeFZkxLFiejzzU7WtnBo+Te6UUhNDW4utMi0/aJO38kPdl+H2Pg0MsYlaV/I2szuZi5mqY7xNAG3tHewtquGTnO6SubK65kH3m5McyYqMeE7JiOOUjHgdIFiNKE3ulFLjhzF2WqyuBM6lNK4yd3jTZYVEQXxGz9K3ztva9m3CaWxpZ9uxSjYerWBzbgXbjlXR0DLw+yo4UFicHsuKjDhOmR7P8ulxOm2X8ipN7pRSY09LQ3cCV3awO4ErPwwtdUM/XmAoJMyChNn2kjgH4mfZJC4iUcd7m8AaWtrYklvJxiMVfHyknB35VbS2D9xeLiosiBXT45ySuXgWp8cQFqx/AtTo0eROKeW/WuqhdL9z2edcZ9tSuOH0RI2Z2jOBS5gFCXPsch06RAH1zW1szq3k4yPlbDxSzs78wTs/TImdZNvKOdWsc5OjCAjQPwTKdzS5U0r5XnMtlB5wErjs7mSu6tjQjxUWYxO2hNmQ6CRyCXNsKVxI+MjHrsa02qZWl2Sugl0F1bQPkszNTo5k5Yx4Tp0Rz4qMeKbEThqlaJXyjCZ3SqnR01TdqxRuH5Tsg5r8oR1HAmybt6RMpwRudndCp9WoagDVja1szqlg41Fbzbq7oJpBcjkyU6JYOTOeVTMTOHVGvA5JovyeJndKqZHX2ggl2XB8j710JnO1hUM7jgTaqtOkTEia131JmK2T1iuPVDe0sinHJnIbj5azt7Bm0GRu3uQoVs1MYNXMeE6dkUC8dn5QY4wmd0qp4TPGzsRwfA8c3wXFuznl6CZ4f4hzowYE2YQtKROS5ncncwmzdDBfNST1zW1syqlg/aEyPjpUTnZxDWaAZE4EFqRGs3JGAitnxnNqRrz2ZFVjniZ3SinPtDbZErjju6F4t70+vhsaK3tsFjHQMQJDbPVpV0lcJiTPt+3hAoO9Gr4an1raOtiRX8WHB8tYf7iMbceqBuwAESCwMC2GlTNsNespGfHEhOt7T40vmtwppXrqHCuuM3nrTOTKDg5hnDixCdvkRZCyyCZwSfNsO7lA/dpRw9fRYcgurmH9oXI+OlzGpqMVA44zFxggLEqLZtVMWzK3IiOe6DBN5tT4pt+ySk1kHe1QdgAKt0Pxru6ErqHc82OERkPKQueyiC2FLSy/8AYIGbAMTymPGGM4VtHAR04yt+FwORX1LQPuM29yFGfMTuSM2bZkLkqTOTXBaHKn1ETR0WHnSy3cBoVb7XXRDmht8PwY8TOdJO4kp1RuoZ3s3qV3am1WliZ26oSU1jaz/nAZ6w+V8+GhMgqqGgfcfmr8JM6YlcjpsxM5fVaC9mZVE54md0qNR8ZAVS4UOElcZyLXXOPZ/iGRXSVx3VWrC3Sye+UVtU2tbDpaYUvnDpWx/3jtgNsnRIRw2qwEWzo3K5FpCTp+oVKuNLlTaqzr7LFa6JLIFW7r09GhX1GpkLYMJi/uTuRip+uMDcpr2jsMuwqqWXeglPcPlLI9r2rAgYPDQwJZOSPeqWpNJDNFZ4BQaiCa3Ck11tSV9CyRK9wG9SWe7RueaBO5KSfb69SlEJ3q3XiVAkpqmlh3sIz3D5Ty4cFSKhta+902OFBYNjWuq93ckqmxBAfqnw2lPKXJnVL+rK0FirbDsY8hb6NN5GoKPNs3LMYmcGlOIpe2DGLSdfYGNSpa2jrYnFvB+wdKWXegjOyi/psEdI41d4bTZu7UGfGEh+jPk1LDpZ8epfxJQwXkbYK8j+HYRlvV2tY0+H4hkbYULm1pd6lc3AxN5NSoyi2v76pqXX+4fMAhSpKiQjl7ThJnz03krDlJOguEUiNIkzulfMUY23s1b2N3yVzpvsH3C5oEqYt7lsolzNY2cmrU1Te38fGRcqd0rpSc8v57XgcHCiumx3P23CQ+NTeJ+alRiP75UMorNLlTarS0tUDxTji2wUnmNnnWVi4uA6augjofEjkAACAASURBVGkrIf1UOxiwDgSsfMAYw77i2q7Suc05lbS09z/N3LT4cD7lJHOrZiUQGarvW6VGg37SlPKWxkqbwHWWyhVsGbyKNSAIUpd0J3NTV0LU5NGJVyk3appa+eBAGe/tL+GDg6Ucr2nud9vwkEBOm5nQVTqXkajjHSrlC5rcKTVSao/D0fch9yPbXq40e/B9QmNg6qlOIrcKpiyHEB2zS/lWTlk9/9pXwr+yj7PpaMWAc7XOmxzVVTq3PCOO0KDAUYxUKeWOJndKDVdzHeSuhyNZcOQ9KNk7+D6x02HaKnuZuspWsWpbOeVjbe0dbMmt5Ll9LdyzJYvDpfX9bhsbHsyZsxP51Nwkzp6bREp02ChGqpTyhCZ3Snmqvc32Xj2SBYffg/xN0NHW//YBQXZg4GmrbPXqtFVaxar8RnVDK1kHSnh3XwlZ+0upbuwcd67v+HMnTYnh0/OSWZ2ZxOL0WAJ1AGGl/Jomd0r1xxgoO+iUzGVBzgcDT98VEGwTuBln2+spy3WOVeVXDpfW8W52Ce9kH2dzbmW/s0KEBQdw5uwk1sxP5tPzkrV0TqkxxmvJnYhMBZ4EUgADPGKM+a2I/AT4GlDqbHqXMeYNZ58fAjcD7cDtxpg3neUXAr8FAoFHjTH3eStuNcF1tpvrTOgGGzA4ZRHMXA0zz4Hpp2kyp/xKa3sHnxyt4F/7bAnd0bL+q1tTY8KYF93GTWuWcdqsBMKCte2cUmOVN0vu2oA7jTFbRSQK2CIibzvrfmOM+ZXrxiKyALgGWAikAe+IyFxn9YPAeUA+8ImIvGqM8aCBk1KD6NFuLgtK9gy8fXS6TeZmnWNL6CKTvR+jUkNQWd9C1oES3skuYd3+Umqb+286sGRqLOfOS+bT85NZkBrN+++/z+p5+p5WaqzzWnJnjCkCipzbtSKSDUwZYJfLgeeMMc3AURE5BJzqrDtkjDkCICLPOdtqcqeGzhg7ndfBt20yl7cJOvqf45LQGJhxVnfpXMIsnfVB+Z1j5Q28uaeYt/YWsyW3kv46t4aHBHLWnETWzEth9bwkkqO0ulWp8WhU2tyJSAawDNgInAHcJiI3AZuxpXuV2MTvY5fd8ulOBvN6LV/p5jFuAW4BSEpKIisra0Sfgxo9dXV1I3r+pKOdmOo9JJZ9TGLZRsKay/rdtkOCqImeR0X8EirjllAXORsTEAgNwO587NtPDWSkz5/qyxhDYZ1h8/E2thxv51ht/wMJJ4QJS5MDWZoUSGZ8ICGBdVBfx94th/v8Q9ZzN7bp+VOdvJ7ciUgk8CLwLWNMjYg8BNyDbYd3D/Br4Csn+jjGmEeARwAyMzPN6tWrT/SQykeysrI44fPX2giH34Xs1+HAP+yAwv1xaTcXMP00YkMiiD2xR5/QRuT8qT6MMezMr+afe4p5c3cxR8rcT/UlAsumxrJmfgpr5ieTmeL5NF967sY2PX+qk1eTOxEJxiZ2zxhjXgIwxhx3Wf8n4HXnbgEw1WX3dGcZAyxXqltjJRx4E7Jfs4ldaz/zXIbFwtwLYc552m5O+bX2DsPmnIquhK6w2v0MJyGBAZw5J5ELF07m0/OTSYwMHeVIlVL+xJu9ZQV4DMg2xtzvsjzVaY8H8Dlgt3P7VWCtiNyP7VAxB9gECDBHRGZgk7prgOu8FbcaY2oKYd/fYd/rkPNh/+PORU+BeRfby/QzIDB4dONUykMtbR2sP1xm29DtOU55fYvb7cJDAjknM5kLFk3mnMwkosL0Pa2UsrxZcncGcCOwS0S2O8vuAq4VkaXYatkc4FYAY8weEXke21GiDfimMaYdQERuA97EDoXyuDFmkC6NalwrPQD7XrNJXcGW/rdLzLTJ3PxLIO1k7Qih/FZjSzvvHyjlzT3FvJN9nNom939SYiYFc+78FC5cNJmz5iTqcCVKKbe82Vv2Q2ypW29vDLDPvcC9bpa/MdB+apzr6IDCbd0JXdmB/redssJJ6C6FxDmjF6NSQ1TT1Mq72SX8c3cxWQdKaGp13ykiMTKUCxbahG7VzASCA3W6OqXUwHSGCuWf2tvsjBD7/m4vtYXutwsIgowzYd4lNqmLThvdOJUagsr6Ft7cU8w/9xTz0aEyWtvdj1kyJXYSFy2azIWLJrNsWpxO96WUGhJN7pR/KT/MjCNPwZavQ22R+22Cw2H2Gph3Kcw9HybFjW6MSg1BXXMbb+8t5tXthXxwsIy2fgahm50cyYULbUK3MC3a4x6uSinVmyZ3yvdaGmDvK7Dtacj9kOnutpkUB5mfsSV0s86B4EmjHaVSHmtqbSdrfwmv7ijkX9klNLe5r3I9aUoMFy6azAULU5idHDXKUSqlxitN7pRvGAOFW2HrU7D7RWiu6btNRBIs/LztEDHtdAjUt6vyX23tHXx0uJxXtxfy1p7ifqf9Wjo1lksWp3LBwslMjQ8f5SiVUhOB/lqq0VVfDjv/akvp3M3jKoGUxZ9M4nnfhjnn65Alyq91dBg251by6o4C3thVTEU/w5bMmxzFpUvSuGxJmiZ0Simv0+ROeV9HOxx5z5bS7X8D2t38AMbPgpNvhCXXsnvLPlbPWz3qYSrlCWMMuwtqeHVHAa/vLKKon4GFpyeEc9mSNC5dksbcFK1yVUqNHk3ulPdU5sC2Z2D7WqhxMydrcDgs+KxN6qad5jIO3b7RjFIpjxwqqeXVHUW8tqOQo2X1brdJiQ7lksW2hG5xeox2ilBK+YQmd2pktTbZ6b+2PQlH17nfZsoKm9At/DyERY9ufEoNQX5lA6/tKOLVHYVkF7lpFwrEhQdz0UmpXLYkjVMz4gnQYUuUUj6myZ0aGUU7bLXrruehqbrv+vAEWHyNTeqS549+fEp5qKK+hdd2FPLK9gK2Hqtyu01ESCAXLJzMpUvTOHN2og4srJTyK5rcqeFrqbdVrlufgOJdfddLAMxaYxO6uRdBUMjox6iUB9raO3j/QCkvbMnnnezjbgcXDgkKYM28ZC5bksY585J16i+llN/S5E4NXUs9fPIofPS/0FDWd31cBiy7AZZcBzFTRj08pTx14HgtL2zJ56WtBZTVNfdZHxggnDUnkcuWpHHeghSiwrT3tlLK/2lypzzXXAef/AnW/w4aynuuCwqD+ZfZUrrpZ0KAVlMp/1Td0MqrOwp4YUs+O/LdNCHAjkX3heXpXHxSKvERWuKslBpbNLlTg2uuhU2PwPrfQ2NFz3UxU+H022HxVTAp1jfxKTWI9g7DuoO22vXtPcdpae87Y0RyVCifO3kKVy5P19kilFJjmiZ3qn9NNTap2/B7aKzsuS5mGpx9p6161bZ0yk8dKqnjhS35vLwtn+M1fatdQwIDOG9BClesSOes2YkEaccIpdQ4oMmd6qupBjY+bJO6pl69BWOnwVnfhSXXalKn/FJ1Yyuv7yzkhS35bOunt+vi9BiuWJ7OZUvSiA3X97FSanzR5E51a6p2kroH+yZ1cRlOUneNTgmm/E57h+GjQ2W8sCWfN/cU09zWt9o1MTKEzy2bwheWpzNvso6vqJQavzS5U9BYBRv/CB//oe8YdXEz4Ozv2TZ1mtQpP3O0rJ4XtuTx0tYCt9OABQcKa+alcMXydD6VmaTj0SmlJgRN7iayxir4+CF7ae6V1MXPtEndSVdBoL5NlP9obmvnn7uLeebjY2zKqXC7zYLUaK5ckc7lS6dob1el1ISjv9oTUWMlbPiDLa1r7jWlUsJsm9QtukKTOuVX8isbWLvxGM9vzqOsrqXP+viIEC5fmsYVy9NZmBbjgwiVUso/6K/3RNJQYateNz7sJqmbA5/6Piz6AgToyPvKP3Q4Q5g8/XEu7+4roaPXxBGBAcI5mclcuSKdczKTCQnSalellNLkbiJoqLA9Xzc+Ai21PdclzoWzvw+LPq9JnfIbFfUt/G1zHs9sPMaxioY+6ydHh3Hdymlcc8pUkqPDfBChUkr5L03uxjNjYMdz8I/v9y2pS8y0JXULP6dJnfILxhi25VXx9IZcXt9VRIubHq9nzUnkhlXTWTMvWcekU0qpfmhyN141VMBrd0D2qz2XJ82zSd2Cz2pSp/xCQ0sbr2wv5KkNuewtqumzPmZSMFcuT+f6VdOZkRjhgwiVUmps0eRuPDr0DvzfN6GuuHtZ/ExY818w/3Kd91X5hUMltTz98TFe3JJPbXNbn/VL0mO4ftV0Ll2cxqQQ/SOilFKe0uRuPGlthLd/DJse7rl8xVfg/J9BiJZ6KN9qbe/g7b3HeWpDLhuOlPdZHxoUwOVL07hh1XQWp+tcxUopNRya3I0XhdvhpVugbH/3sogkuPxBmHuB7+JSCiiqbuTZTXk8t+kYJbV953idkRjB9SunccXydJ0OTCmlTpAmd2NdRzt89AC893PocKnayrwYLvtfiEj0XWxqQjPG8PGRCv6y/ijvZJfQ3msck8AA4bz5Kdywajqnz0ogIEB8FKlSSo0vmtyNZZU58PLX4diG7mXBEXDRfbDsRhD9sVSjr73D8NqOQh5Zd4RdBdV91idHhXLNqdO49tSppMZM8kGESik1vmlyNxYZA9vXwj/+o+e4demnwucftp0nlBpl9c1tPL85jwc/aKSscVuf9afNTODG06Zz3oIUneNVKaW8SJO7saa+HF6/A7Jf614mgbD6h3Dmt3XKMDXqSmqbeHJ9Lk99nEt1Y2uPdaFBAVyxPJ0vn5HB7OQoH0WolFITi2YCY8nBd+CVb0Dd8e5lCbPh84/AlOW+i0tNSIdK6nj0gyO8tLWAlvaeAw7HhQdz02kZ3HjadBIjQ30UoVJKTUya3I0FLQ3w9n/BJ3/quXzFzXD+PTrEiRo1xhg+yankkXWHeSe7pM/66QnhnJ3Sxl3XfFrHplNKKR/R5M7fFW5zhjg50L0sItkZ4uR838WlJpT2DsObe4p5ZN0RtudV9Vm/dGost549k/MXTuaDde9rYqeUUj6kyZ2/6miHD++HrPt0iBPlM40t7bywJY9HPzxKbnlDn/Xnzk/hlrNnckpGHKK9s5VSyi9ocuePKo7Cy7dC3sbuZTrEiRpF5XXNPLEhl6c25FDZ0LOTREhgAJ8/eQpfPWsms5MjfROgUkqpfmly50+Mge3POEOc1HUv1yFO1Cg5WlbPox8c4YUt+TS39ewkETMpmBtWTeOLp2eQHBXmowiVUkoNRpM7f1FfDq/dDvte714WEASf+oEOcaK8bkuu7STx1t7jmJ4TSZAeN4mbz5zBVSumEhGq70OllPJ3+k3tD6rz4bELoCa/e1nCHGeIk5N9F5ca9zbnVHD/2wdYf7i8z7pFU6K55exZfGbRZIJ00GGllBozNLnzteY6ePaanondKV+F8+6BkHDfxaXGte15Vdz/9gHWHSjts251ZhK3nD2T02YmaCcJpZQagzS586WODttxoniXvR8QBFc9CfMu9m1catzaXVDNb94+wL/29RyjLjBAuHxpGreePYvMyTqThFJKjWWa3PnSuz/t2cbu4vs1sVNekV1UwwPvHODNPcd7LA8Q+OyyKdz+6TlkJOpg2EopNR5ocucr29fCh7/pvn/abbD8i76LR41LB4/X8sC/DvL3nUU9lovAJYvTuGPNHB3ORCmlxhlN7nwhdwO8env3/bkXwnk/9V08atw5WlbPb985wCs7Cvv0fr1o0WS+de5crX5VSqlxSpO70VZxFP56PXQ4A8MmL4QvPAoBOl2TOnHHyhv433cP8vK2Ato7emZ1585P4dvnzWFhWoyPolNKKTUaNLkbTU3VtmdsgzPsREQSXPcchGoJijox+ZUNPPjeIf62OZ+2Xknd6swkvnPeXBanx/ooOqWUUqPJa8mdiEwFngRSAAM8Yoz5rYjEA38FMoAc4CpjTKXYMRd+C3wGaAC+ZIzZ6hzri8CPnEP/zBjzhLfi9pr2NnjhK1C6z94PDIFr1kLsNN/Gpca04uomHnzvEM99cozW9p5J3VlzEvnWuXNZPj3OR9EppZTyBW+W3LUBdxpjtopIFLBFRN4GvgT8yxhzn4j8APgB8B/ARcAc57ISeAhY6SSDPwZWYJPELSLyqjGm0ouxj7y3/hMOvdN9//IHYeqpvotHjWkltU08lHWYZzYeo6XXNGErZ8TznfPmsnJmgo+iU0op5UteS+6MMUVAkXO7VkSygSnA5cBqZ7MngCxscnc58KQxxgAfi0isiKQ6275tjKkAcBLEC4FnvRX7iPvkMdj4x+77Z30XFl/lu3jUmFVe18zD647w5IYcmlp7JnXLp8dx53lzOW2WDj6slFIT2ai0uRORDGAZsBFIcRI/gGJstS3YxC/PZbd8Z1l/y3s/xi3ALQBJSUlkZWWNWPwnIq5iO4t33k3nT21J0unsDTgD/CQ+f1RXV+c3589fNLQa3jjaytu5rTS391w3IyaAz88OZlFiMy35u3k/3/0xRouev7FLz93YpudPdfJ6cicikcCLwLeMMTWuJQrGGCMipt+dh8AY8wjwCEBmZqZZvXr1SBz2xJQdhEdvApwSltSlJH/5RZJ1WrEBZWVl4Rfnzw+0dxhe2JLHL9/cT1lda491C9Oi+c55c/n0vGS/KqnT8zd26bkb2/T8qU5eTe5EJBib2D1jjHnJWXxcRFKNMUVOtWvnPEgFwFSX3dOdZQV0V+N2Ls/yZtwjoqEC1l5te8gCRKXCtc/qfLHKY5uOVnD3a3vYU1jTY3lmShTfPm8uFyxM8aukTimllH/wZm9ZAR4Dso0x97usehX4InCfc/2Ky/LbROQ5bIeKaicBfBP4uYh0dvk7H/iht+IeEe2t8PxNUHHY3g+aZBO76DTfxqXGhIKqRv77jWxe7zWrRGpMGD+4aB6XLk4jIECTOqWUUu55s+TuDOBGYJeIbHeW3YVN6p4XkZuBXKCzZ8Eb2GFQDmGHQvkygDGmQkTuAT5xtvtpZ+cKv2QM/P1OyPmge9nnH4G0Zb6LSY0JjS3t/PH9w/zx/cM0u/SADQ0K4OufmsXXPzWLSSE62LVSSqmBebO37IdAf8ULa9xsb4Bv9nOsx4HHRy46L/r4D7DVZRi+T/8/WHCZ7+JRfs8Yw2s7i/jvN7Ipqm7qse6Sxan88DPzmRI7yUfRKaWUGmt0hoqRdOBNePM/u+8vvhrOutN38Si/tyu/mrtf28Pm3J7DNi5Mi+bHly7k1BnxPopMKaXUWKXJ3Ug5vsfOQIHT+XfqSrj0f0EbvCs3Smqb+NWb+/nblnyMS3/xxMgQvnt+JleumEqgtqtTSik1DJrcjYS6Ulh7DbTU2fsx0+DqZyA4zLdxKb/T3NbOXz7K4XfvHqKuua1reXCg8OUzZnDbp2cTHRbswwiVUkqNdZrcnajWJnjuOqg+Zu+HRMJ1z0Fkkm/jUn7FGMO/skv42d/3klPe0GPdmnnJ/OfF85mZFOmj6JRSSo0nmtydCGPgtdshf5O9LwFwxeOQstC3cSm/cvB4LT99fS8fHCzrsXxWUgT/75IFrM5M9lFkSimlxiNN7k7EB7+GnX/tvn/+z2DuBb6LR/mVqoYWHnjnIE99nEt7R3fDuuiwIL593lxuWDWd4MAAH0aolFJqPNLkbrj2vgLv3tN9/+Qvwqpv+C4e5Tfa2jt4dtMxfv32AaoauqcMCxC4buU0vnNeJvERIT6MUCml1Himyd1wFG6Dl27tvp9xFlz8a+0Zq1h/uIy7X93L/uO1PZavmhnPjy9dyPzUaB9FppRSaqLQ5G6oagrh2WuhrdHej58FVz0JgdrDcSKraWrl53/P5rlP8nosT4+bxI8uns8FCyfrPLBKKaVGhSZ3Q9HSYBO7WmfOz7AYuO6vEK4DzU5k7x8o5Qcv7uwxu0R4SCDfPGc2N585g7BgnTJMKaXU6NHkbije+B4UOdPkSiBc+QQkzvFtTMpnaptauddNad2FCyfzk8sWMjlGxzlUSik1+jS581RrE+x4tvv+Z34Js87xXTzKp9Y5pXWFLqV1ceHB/PTyRVyyOFWrYJVSSvmMJneeqjoGpt3ejpkGp9zs23iUT/RXWnfRosnc89lFJEaG+igypZRSytLkzlNVud2346b7Lg7lM1pap5RSaiwYNLkTkdOAG4CzgFSgEdgN/B142hhT7dUI/UVlTvftuAxfRaF8QEvrlFJKjSUDJnci8g+gEHgFuBcoAcKAucA5wCsicr8x5lVvB+pzPZI7LbmbKLS0Timl1FgzWMndjcaYsl7L6oCtzuXXIpLolcj8TY/kbobPwlCjQ0vrlFJKjVUDJnduErthbTMuuLa5i9WSu/FMS+uUUkqNZR51qBCRVcDvgPlACBAI1BtjJsZcSsZApXaoGO9qm1r5+RvZPLup77h193x2EUlRWlqnlFLK/3naW/b3wDXA34AVwE3YdncTQ2MlNNfY28HhEJHk23jUiNPSOqWUUuOFx0OhGGMOiUigMaYd+LOIbAN+6L3Q/EjvKln9oR83tLROKaXUeONpctcgIiHAdhH5H6AICPBeWH5Gh0EZl7S0Timl1HjkaXJ3I7ad3W3At4GpwBe8FZTf0fZ244qW1imllBrPPErujDGd2U0jcLf3wvFTWnI3buwprObrT28hr6Kxa5mW1imllBpPBhvE+HljzFUisgswvdcbYxZ7LTJ/osOgjAsvb8vnBy/uormto2uZltYppZQabwYrubvDub7E24H4NS25G9Na2jr4+RvZ/GV9TteyyNAg7v3cIi5bkqaldUoppcaVwQYxLnJulgGNxpgOEZkLzAP+4e3g/EJHO1S5tM2Knea7WNSQldQ08Y1ntrI5t7Jr2aykCB6+cQWzkyN9GJlSSinlHZ52qFgHnCUiccBbwCfA1cD13grMb9QUQkervR2RBKGaEIwVm3Mq+LdntlJa29y17KJFk/nllUuIDPV4FCCllFJqTPH0F06MMQ0icjPwB2PM/4jIdm8G5je0vd2YY4zhqY9z+elre2nrsE1FAwS+f+E8bj17plbDKqWUGtc8Tu5E5DRsSd3NzrJA74TkZ3q0t9Pkzt81tbZz18u7eGlrQdeyuPBgfnftyZw5J9GHkSmllFKjw9Pk7g7sbBQvG2P2iMhM4D3vheVHeoxxl+GzMNTg8ioauPWpLewtquladtKUGB664WTS48J9GJlSSik1egZN7kRkKZAM/IcxJhvAGHMEuN3LsfkHrZYdE94/UMrtz26jurG1a9mVy9O557OLCAueGIXMSimlFAwyhZiI/BfwPHY2ir+LyNdGJSp/osOg+LWODsOD7x3iS3/e1JXYBQcK935uEf9zxWJN7JRSSk04g5XcXQ0sdTpTJAD/BP7k/bD8iE495rdqm1q58/kdvLX3eNeylOhQHrphOSdPi/NhZEoppZTvDJbcNRtjGgCMMeUiMmBJ37jT2gh1xfa2BEJ0um/jUV0OHq/l1qe2cKSsvmvZqTPiefC6k3W2CaWUUhPaYMndTBF51bktwCyX+xhjLvNaZP6g6lj37Zh0CNSx0fzBG7uK+O7fdtDQ0t617OYzZ/CDi+YRHDix/n8opZRSvQ2WrVze6/6vvBWIX9L2dn6lrb2DX761n4ffP9K1LCw4gF98YTGXL53iw8iUUkop/zHY9GPvj1Ygfknb2/mN8rpmbn9uGx8dKu9aNj0hnD/esJz5qdE+jEwppZTyLwMmdyLyGvAI8E9jTGuvdTOBLwE5xpjHvRahL7mW3OkwKD6zM7+Kf3t6KwVVjV3LzslM4oGrlxETHuzDyJRSSin/M1i17NeA7wAPiEgFUAqEARnAYeD3xphXvBqhL1XpAMa+9vwnefzold20tHV0LbtjzRzuWDOHgACdRkwppZTqbbBq2WLg+8D3RSQDSAUagQOdvWjHNW1z5zPNbe3c/dpe1m7s7tQSFRbEA1cvZc38FB9GppRSSvk3j7t/GmNygByvReJvjNGpx3yksc1w02Ob2Hi0omtZZkoUD9+4nIzECB9GppRSSvk/HdujP42V0FJrbwdHQHiCb+OZIKobW/nVJ00cru4uGL50SRq/+MJJhIfo21UppZQajP5a9qfyaPftuAwQbd/lbZX1Ldz4+EYOV3e3r/veBZl8Y/UsRF9/pZRSyiMeJ3ciMgmYZozZ78V4/IcOgzKqSmubueHRjew/Xtu17O7LFvLF0zN8F5RSSik1Bnk0nL+IXApsx84ti4gsdZ2pop99HheREhHZ7bLsJyJSICLbnctnXNb9UEQOich+EbnAZfmFzrJDIvKDoT7BYdPOFKOmuLqJqx/Z0JXYCfCLL5ykiZ1SSik1DJ7O1fQT4FSgCsAYsx2YMcg+fwEudLP8N8aYpc7lDQARWQBcAyx09vmDiASKSCDwIHARsAC41tnW+1yHQdEx7rwmv7KBqx7ewJFSO0dsYIDwtcWhXH3KNB9HppRSSo1NnlbLthpjqnu1ezID7WCMWecMn+KJy4HnjDHNwFEROYRNJgEOGWOOAIjIc862ez087vBpyZ3X5ZTVc92fPqawugmAoADhd9cuY1L5xKj5V0oppbzB05K7PSJyHRAoInNE5HfA+mE+5m0istOpto1zlk0B8ly2yXeW9bfc+7TNnVcdPF7LVQ9v6ErsQgID+OMNy7nopFQfR6aUUkqNbZ6W3P078J9AM7AWeBP42TAe7yHgHmyp3z3Ar4GvDOM4fYjILcAtAElJSWRlZQ3/YKads6uOdWW+63bl0hF4/IRjVNaxmnZ+ubmJ2hZ7PyQAbl8WQlBJNlkl2dTV1Z3Y+VM+pedv7NJzN7bp+VOdBk3unHZvfzfGnINN8IbNGNOVIYnIn4DXnbsFwFSXTdOdZQywvPexH8HOg0tmZqZZvXr18AOtOgbvt9vbEcmcvcZd00E1HDvyqrjj8U1diV1ESCCPfekUVs3sHkcwKyuLEzp/yqf0/I1deu7GNj1/qtOg1bLGmHagQ0RiTvTBRMS1zu1zQGdP2leBa0QkVERmAHOATcAnwBwRmSEiIdhOFwP20h0RPdrbePdW6AAAIABJREFUaZXsSNmcU8ENj26kurEVsNOJPXnzyh6JnVJKKaVOjKfVsnXALhF5G6jvXGiMub2/HUTkWWA1kCgi+cCPgdUishRbLZsD3OocZ4+IPI/tKNEGfNNJKhGR27DVwIHA48aYPUN5gsOi046NuPWHy/jqE5tpaLElorHhwTx980oWTTnh/wxKKaWUcuFpcveSc/GYMeZaN4sfG2D7e4F73Sx/A3hjKI99wnQYlBGVtb+EW5/aQnObnXkiMTKEp7+6knmTo30cmVJKKTX+eJTcGWOe8HYgfkWHQRkxb+0p5ra122hpt4ldSnQoz3x1FbOTI30cmVJKKTU+eZTcichR3IxrZ4yZOeIR+QMdBmVEvL6zkG89t522DvvWmRI7ibVfW8n0hAgfR6aUUkqNX55Wy65wuR0GXAnEj3w4fkJL7k7Yi1vy+d4LO3DyOqYnhLP2a6uYEjvJt4EppZRS45xHgxgbY8pdLgXGmAeAi70cm2+0NEB9ib0dEATRozNm8niyduMx7vxbd2I3KymC5289TRM7pZRSahR4Wi17ssvdAGxJnqelfmOLa2eKmKkQEOi7WMagP390lLtf654dbt7kKJ7+6koSI0N9GJVSSik1cXiaoP3a5XYbdhiTq0Y8Gn+g7e2G7aGsw/zin/u67i9Oj+HJr5xKbHiID6NSSimlJhZPe8ue4+1A/IZrezsdBsUjxhgeeOcgv/3Xwa5ly6fH8ecvn0J0WLAPI1NKKaUmHo/a3InIHSISLdajIrJVRM73dnA+UaUDGA+FMYb7/rmvR2J32swEnvzKqZrYKaWUUj7gUXIHfMUYUwOcDyQANwL3eS0qX9KpxzzW0WG4+7W9PPz+ka5lZ89N4s9fPoWI0PHZJFMppZTyd54md+JcfwZ40pkCTAbYfuzSqcc89os39/GX9Tld989bkMKfblpOWLB2QlFKKaV8xdPkbouIvIVN7t4UkSigw3th+YgxvaYey/BZKP7ule0FPUrsLl6cyh+uP5nQIE3slFJKKV/ytO7sZmApcMQY0yAi8cCXvReWjzSUQ0udvR0SBeHjd5zmE7G7oJr/eHFn1/0185L57dVLCQr09L+CUkoppbzF01/j04D9xpgqEbkB+BFQ7b2wfKT3MCgyPmueT0R5XTO3PrWFplZbcDszKYLfXKOJnVJKKeUvPP1FfghoEJElwJ3AYeBJr0XlK5VHu29re7s+Wts7uG3tNgqqGgGICg3iTzet0F6xSimllB/xNLlrM8YY4HLg98aYB4Eo74XlIz3a22lP2d5+/kY2G46UA7ZQ84FrljIrKdLHUSmllFLKladt7mpF5IfYIVDOEpEAYPwV1+gwKP16YUs+f/4op+v+d86dy5r5Kb4LSCmllFJueVpydzXQjB3vrhhIB37ptah8RYdBcWtHXhV3vbyr6/6FCyfzzXNm+zAipZRSSvXHo+TOSeheBDpnfy8DXvZWUD6jU4/1UVprO1C0tNkOFHNTIvnVVUsICNDOJkoppZQ/8nT6sa8BLwAPO4umAP/nraB8or0NqvO778dO810sfqKlrYNvPLOF4pomAKLDbAeKSJ19QimllPJbnlbLfhM4A6gBMMYcBJK9FZRP1OSDabe3I1MgJNy38fiBn/7/9u48vorq/v/465OETVkSllIkaEAoWzZWWQwGKagUQS0Wv7JpC35dU+lPCq2VunxrsVJRi8tXqwIuqMUvSAF3SMG6IAgVkB0iBJA9QMCwnt8fd3K5ISshyU3ufT8fjzwyM3dmzmdmuObjOXPOmbuarzIOABBh8LebO3JJgwuDHJWIiIgUpaTJ3THn3PHcFTOLAlz5hBQket8ujxlLtvLaF1v96+OubsMVP2kUxIhERESkJEqa3P3LzH4P1DKzvsA/gH+WX1hBoGFQ/JZ9t58J767yr1+bdBG39WoRxIhERESkpEqa3I0D9gArgf8G5uObpSJ05BkGJS5YUQTdrkM53P7a15w45auYbdekLn/5eSKm2TpERESqhGLfjDezSGC1c64N8GL5hxQkZ089FoaOnTzFf7+6jD2HjwEQc0E1/nd4J2pVjwxyZCIiIlJSxdbcOedOAevMLLS7j4Z5zZ1zjgdmr2LFtiwAIiOMZ27uSLP66lgiIiJSlZR0TIsYYLWZLQGO5G50zg0sl6iCIczfuXvti+94e+mZoWDu79+WHi0bBjEiERERKY2SJncPlGsUwXYsG47s8S1HVIO6FwU3ngr25eZ9PPTPb/3rN3Rsyq0944IXkIiIiJRakcmdmdUEbgda4utM8ZJz7mRFBFahss4M+UF0M4gIn3fMtmf9wJ2vf83J074OFImx9Xj0+gR1oBAREamiinvnbhrQGV9idw3w13KPKBjCdNqxnBOn+O9Xl7LviG8Iw4a1q/P8sE7UrBY+ya2IiEioKa5Ztp1zLgHAzF4ClpR/SEGQFX4DGDvn+N3/rWTV9kMAREUYzw7txEXRtYIcmYiIiJyP4mruTuQuhGRzbK48PWXDo+bupU+3MGv5dv/6Hwe2p2vz+kGMSERERMpCcTV3SWZ2yFs2fDNUHPKWnXOubrlGV1HCbOqxTzfs5dH5a/zrN3VpxrDLQnukGxERkXBRZHLnnAuPl6/CaBiUbfuPcveMr/H6T9Dh4mgeGtReHShERERCREmnHwtdzoXNAMZHj59k9PSlZB31tbb/qE4Nnh/WiRpR4ZHDi4iIhAMld0f2womjvuUadaFWTHDjKSfOOcbO/Ia13x8GoHpkBM8P70TjujWDHJmIiIiUJSV3Z3emCNHmyef/tZl53+z0rz88qD0dLw7NRFZERCScKbkLg/ft0tft5i8frPWvD+92CTd1VQcKERGRUKTk7sCWM8sh+L7dtv1HSZuxHOd1oOgaV58HBrQLblAiIiJSbpTchfgwKI/M/ZZDOb4hCpvUq8kzQztSPUqPXUREJFTpr3wITz322ca9fPjtLv/6lJs70qhOjSBGJCIiIuVNyV2ITj126rTj4bnf+tdv6NiUTpeoA4WIiEioC+/k7tQJOJh5Zj06dDoZvPnVVv+wJ7WqRTLu6jZBjkhEREQqQngndwczwZ32LddpAtVCY8y3Qzkn+OuH6/3rd6ZeqvHsREREwkR4J3chOgzKlAUb2X/kOABNo2sxuleLIEckIiIiFSW8k7sQnHZsy94jvPLvM8O7jL+mDTWraXoxERGRcBHmyV1gZ4rQqLl7dP4aTpzyDWrX+ZIYBiQ2CXJEIiIiUpHKLbkzs5fNbLeZrQrYVt/MPjKzDd7vGG+7mdnTZrbRzL4xs44Bx4z09t9gZiPLNMgQGwbl3xv38lHA0CcTrm2Hheh0aiIiIlKw8qy5mwpcfda28cAnzrlWwCfeOsA1QCvv5zbgOfAlg8AfgcuArsAfcxPCMhFCw6CcPHWah/95ZuiTn3eMJTE2OogRiYiISDCUW3LnnFsE7D9r8yBgmrc8DbguYPt05/MFEG1mTYCrgI+cc/udcweAj8ifMJZennfuqnbN3ZtfbWPdLt/QJxdUj+S3V7cOckQiIiISDFEVXF5j59xOb/l7oLG33BTYFrBfpretsO35mNlt+Gr9aNSoEenp6UUGEnnyKClH9wFw2qJY9PU6sI3nci2VxpETjscWHfWvX31JBGu+/oI1QYzpfGRnZxf7/KTy0vOruvTsqjY9P8lV0cmdn3POmZkrw/O9ALwA0Lp1a5eamlr0Ad+vgk99ixExl5Dau09ZhVLh/jTvWw6f8PWQbRpdi0dHXFGle8imp6dT7POTSkvPr+rSs6va9PwkV0X3lt3lNbfi/d7tbd8ONAvYL9bbVtj28xci79tt2XuEqZ9l+Nd/119Dn4iIiISzik7u5gC5PV5HAu8GbB/h9ZrtBhz0mm8/APqZWYzXkaKft+38hcj7dn+ad2boky5xMfwsQUOfiIiIhLNya5Y1sxlAKtDQzDLx9XqdCLxtZr8CvgN+4e0+H+gPbASOArcCOOf2m9kjwFfefg87587upFE6B6p+zd2nG/by8Rrf0CdmMGFAew19IiIiEubKLblzzv1XIR/le7nNOeeAuwo5z8vAy2UYmk8Vn3rs5KnTPDI379AnCbH1ghiRiIiIVAbhO0NFFZ96bMZZQ5+MvUpDn4iIiEi4JnfOVempxw7+cIInPlznX7+rd0sa160ZxIhERESksgjP5C57N5z8wbdcox7UKrtJLyrC059s4MDRE4Bv6JNfXd48yBGJiIhIZRGeyV1W1a2127wnm2kBQ5/8vn9bDX0iIiIifuGZ3FXhYVAenb+Gk6fPDH3SP+HHQY5IREREKpMwTe6q5jAoizfs4eM1vnGfNfSJiIiIFCRMk7uMM8tVZBiUs4c+GayhT0RERKQA4Znc5Xnnrmp0RpixZCvrd2UDcKGGPhEREZFChGdyV8WGQTl49ARPfLTev35n75b8SEOfiIiISAHCL7k7dQIOZXorBvWaBTWcknhKQ5+IiIhICYVfcndwG7jTvuU6TaBa5a4B27Qnm+mfZ/jXNfSJiIiIFCX8krsqNu3Yo/PODH3SNa6+hj4RERGRIoVhcld13rdbtH4Pn6wNGPrk2nYa+kRERESKFIbJXcaZ5Uo8DMrZQ5/c2CmW+KYa+kRERESKFn7JXVbVGMD4jSVb2bD7zNAn9/XT0CciIiJSvPBL7qrA1GMa+kRERERKKwyTu8pfc/fkJ+vJ8oY+iY3R0CciIiJScuGV3OUcgh/2+5Yja0DtytfzdOPubF79/EwCqqFPRERE5FyEV3IX+L5d9MUQUfku/0/zvj0z9Enz+lwTX/kSUBEREam8Kl92U54q+TAo6et2s3DdHsAb+mSAhj4RERGRcxNmyV3GmeVK9r7dyVOn+Z95a/zrGvpERERESiO8krs8zbKVq+bu9S+3sjFw6JOrNPSJiIiInLvwSu4q6TAoh3NOMPnjM0Of3HVlS35UR0OfiIiIyLkLs+Sucg6D8v6q7/1DnzSrX4tf9tTQJyIiIlI64ZPcOVdpm2XnrdzpXx7e7RINfSIiIiKlFj7JXfYuOJnjW64ZDbWigxuPJ+vocT7dsNe/fk18kyBGIyIiIlVd+CR3lfR9uw+/3eUf1y6pWTTN6l8Q5IhERESkKguj5K5yvm8375szTbIDElRrJyIiIucnjJK7jDPLleR9uwNHjvPvjQFNsgmajUJERETOT/gkd1mVr+buw2+/9zfJdrg4mtgYNcmKiIjI+Qmf5K4STj02N6BJ9mdqkhUREZEyEEbJXcaZ5ZjgjyO3/8hxPtu0z7/eX8mdiIiIlIHwSO5OHodD270Vg3qxQQ0H4IPV33PKa5LtdEkMF0XXCnJEIiIiEgrCI7k7uA3wJVLUvQiiagQ1HMjbS1ZNsiIiIlJWwiO5y9MkGxesKPz2ZR/js01nesmqSVZERETKSvgld5VgGJT3V3+P1yJLl7gYflyvZnADEhERkZARHsldJRsGRU2yIiIiUl7CI7mrRFOP7Tl8jC82+3rJmsE1Su5ERESkDIVJcld5au7yNsnWp3FdNcmKiIhI2QmT5C7jzHKQ37mb980O//KARNXaiYiISNkK/eQu5yDkZPmWo2pC7cZBC2X34Ry+3LIf8DXJXh2vuWRFRESkbIV+chfYJBt9MUQE75LfX/U9zmuSvax5fX5UR02yIiIiUrbCILnLOLMc5Pft8swlm3hRECMRERGRUBX6yV3gMChBfN9u16EcvsrwNclGGFzdXk2yIiIiUvZCP7mrJMOgvLdyp79JtluLBjSqE/wp0ERERCT0BCW5M7MMM1tpZivMbKm3rb6ZfWRmG7zfMd52M7OnzWyjmX1jZh3PqbBKMgzKvJWBTbLqJSsiIiLlI5g1d72dc8nOuc7e+njgE+dcK+ATbx3gGqCV93Mb8Nw5lVIJhkH5/mAOX2UcANQkKyIiIuWrMjXLDgKmecvTgOsCtk93Pl8A0WZWsqqv06cha+uZ9SA1y84PqLXrcWlDGtRWk6yIiIiUj6ggleuAD83MAf/rnHsBaOycy82CvgdyB6RrCmwLODbT27YzYBtmdhu+mj0aNWpEeno61Y/to8epYwCciKrDv79YXl7XU6Q3vvjBv9yyxiHS09ODEkdVkZ2drXtUhen5VV16dlWbnp/kClZyd7lzbruZ/Qj4yMzWBn7onHNe4ldiXoL4AkDr1q1damoqfPc5fO77vNqPWpKamloWsZ+THVk/sPH9BQBERhhpN1xB/QurV3gcVUl6enpQnpWUDT2/qkvPrmrT85NcQWmWdc5t937vBmYBXYFduc2t3u/d3u7bgWYBh8d624pXCd63y9sk20CJnYiIiJSrCk/uzOxCM6uTuwz0A1YBc4CR3m4jgXe95TnACK/XbDfgYEDzbdGygt9TNrCXrOaSFRERkfIWjGbZxsAsM8st/w3n3Ptm9hXwtpn9CvgO+IW3/3ygP7AROArcWuKS8gyDUvE1d5kHjrJ8q29e26gIo1879ZIVERGR8lXhyZ1zbjOQVMD2fUCfArY74K5SFRbkqcfeW/m9f7lny4bEqElWREREylllGgql7AV56rG5GrhYREREKljoJncnj8GhHd6KQb1mRe5e1rbtP8p/tvmaZKtFGlepSVZEREQqQOgmd1nb8A2nB9SLhaiKbRIN7CV7ecuG1LugWoWWLyIiIuEpdJO7IA+DMj9Pk+xFFV6+iIiIhKfQTe6yMs4sV3Bnim37j/KfzIOAr0m2b7vGxRwhIiIiUjZCN7nL01O2YmvuAse269WqEfVqqUlWREREKkYIJ3fBG8B43jfqJSsiIiLBEcLJXcaZ5Qp85+67fUdYud3XJFs9MoKfqklWREREKlDoJndBmnosT5PsTxpRt6aaZEVERKTihGRyZ+405Phqz4iqBbV/VGFlBzbJai5ZERERqWghmdxFuBNnVqIvBt88tuVuy94jrN5xCIDqURH0aVtxSaWIiIgIBGFu2YpgpwOSuwpskg0c2y71J42ooyZZESnEiRMnyMzMJCcnJ9ih+NWrV481a9YEOwwpJT2/yq9mzZrExsZSrVr55gchmdxFnD55ZqUCh0GZq16yIlJCmZmZ1KlTh7i4OKyCWheKc/jwYerUqRPsMKSU9PwqN+cc+/btIzMzk+bNm5drWaHZLBuEmrtNe7JZszOwSVa9ZEWkcDk5OTRo0KDSJHYiUr7MjAYNGlRIbX1IJneW5527iqm5mx9Qa9e7dSNq1wjJSlERKUNK7ETCS0V950MyucvbLBtXIWXO01yyIhLCRo0axbfffgvAo48+6t+ekZFBfHx8mZSxYsUK5s+fX+L9s7KyePbZZ/3rO3bsYPDgwWUSS3p6OgMGDChyn7PjnTNnDhMnTiyT8kuif//+ZGVlAVC7du0SH1eW96koPXr0AHz/Rt54441yLWvKlCm0bNkSM2Pv3r3+7c450tLSaNmyJYmJiXz99df+z6ZNm0arVq1o1aoV06ZN829ftmwZCQkJtGzZkrS0NJxz5Rp7eQjR5C6wWbb8a+427j7M2u8PA1AjKoI+bdRLVkRCy9///nfatWsH5E3uytL5JncXXXQRM2fOLI/QCnR2vAMHDmT8+PEVVv78+fOJjo4+5+Mq6j599tlnQMUkdz179uTjjz/mkkvy/s1/77332LBhAxs2bOCFF17gjjvuAGD//v089NBDfPnllyxZsoSHHnqIAwcOAHDHHXfw4osv+o97//33yzX28hCSyR14WfYFDaBG+b9cOu+b7/3LV7b5EReqSVZEqpB//OMf/OY3vwHgqaeeokWLFgBs3ryZnj17ApCamsrSpUsZP348P/zwA8nJyQwdOhSAU6dOMXr0aNq3b0+/fv344YcfAF/y061bNxITE7n++uv9fzxzzwWwd+9e4uLiOH78OBMmTOCtt94iOTmZt956K0+Mq1evpmvXriQnJ5OYmMiGDRsYP348mzZtIjk5mbFjx+apRZw6dSrXXXcdffv2JS4ujilTpvDEE0/QoUMHunXrxv79+wuN5WxLliyhe/fudOjQgR49erBu3boC4506dSp333034EtorrzyShITE+nTpw9bt24F4JZbbiEtLY0ePXrQokWLEiVZM2bMICEhgfj4eMaNG+ffHhcXl6eW6mzjx4/nmWee8a8/+OCDTJo0Kc99Kui+AkyfPp3ExESSkpIYPnx4kde0a9curr/+epKSkkhKSvIndbm1iePHj2fx4sUkJyczefJkevXqxYoVK/xxXX755fznP/8p9j4UpUOHDgU+u3fffZcRI0ZgZnTr1o2srCx27tzJBx98QN++falfvz4xMTH07duX999/n507d3Lo0CG6deuGmTFixAhmz56d77wFXXNGRgZt2rRh6NChtG3blsGDB3P06NHzuq7SCu0spILet5u3cod/Wb1kReRcxY2fV27nzpj4s2L3SUlJ4S9/+QsAixcvpkGDBmzfvp3FixfTq1evPPtOnDiRKVOm+P84Z2RksGHDBmbMmMGLL77IL37xC9555x2GDRvGiBEj+Nvf/sYVV1zBhAkTeOihh3jyyScLjKF69eo8/PDDLF26lClTpuT7/Pnnn+fXv/41Q4cO5fjx45w6dYqJEyeyatWqPLEEWrVqFcuXLycnJ4eWLVvy2GOPsXz5csaMGcP06dO59957i703AG3atGHx4sVERUXx8ccf8/vf/5533nknX7xTp071H3PPPfcwcuRIRo4cycsvv0xaWpo/Sdi5cyeffvopa9euZeDAgUU2ke7YsYNx48axbNkyYmJi6NevH7Nnz+a6664rNu4hQ4Zw7733ctdddwHw9ttv88EHH3Dq1Kki7+vq1av5n//5Hz777DMaNmzoT4QLu6a0tDSuuOIKZs2axalTp8jOzs4Tx8SJE5k0aRJz584FoH79+kydOpUnn3yS9evXk5OTQ1JSUp5j1q1bx5AhQwq8rvT09BLXWG7fvp1mzZr512NjY9m+fXuR22NjY/NtP1tB13zgwAHWrVvHSy+9RM+ePfnlL3/Js88+y3333VeiWMtSiNbceSrgfbv1uw6zfpfvH3LNahFcqSZZEalifvzjH5Odnc3hw4fZtm0bN998M4sWLWLx4sWkpKQUe3zz5s1JTk4GoFOnTmRkZHDw4EGysrK44oorABg5ciSLFi0qdYzdu3fn0Ucf5bHHHuO7776jVq1axR7Tu3dv6tSpQ6NGjahXrx7XXnstAAkJCfkSwaIcPHiQG2+8kfj4eMaMGcPq1auLPebzzz/n5ptvBmD48OF8+umn/s+uu+46IiIiaNeuHbt27SryPF999RWpqak0atSIqKgohg4dWuL72KFDB3bv3s2OHTv4z3/+Q0xMTJ6EBgq+rwsWLODGG2+kYcOGgC8ZK+qaFixY4G/ujIyMpF69ekXGdeONNzJ37lxOnDjByy+/zC233JJvn9atW7NixYoCf0rTFF3WCrvmZs2a+Wu7hw0blue5V6QQT+7Kv+YucLqxPm0ac0H10K4MFZHQ1KNHD1577TVat25NSkoKixcv5vPPP/f/oSpKjRo1/MuRkZGcPHmyiL0hKiqK06dPA5R4WIibb76ZOXPmUKtWLfr378+CBQvOKa6IiAj/ekREhD/GksTywAMP0Lt3b1atWsU///nP8x7KIjCu8n5Z/8Ybb2TmzJm89dZbBdaElea+nq8LLriAvn378u677/L222/7m/cDrVu3juTk5AJ/cjuRlETTpk3Ztm2bfz0zM5OmTZsWuT0zMzPf9pI6uzdssHrEh3YmUs7Nss65s3rJqklWRM5dSZpOy1tKSgoPPPAAf/zjH+nQoQMLFy6kVq1aBdbCVKtWjRMnThQ5yn69evWIiYnx1/69+uqr/lq8uLg4li1bRteuXfO8c1anTh0OHz5c4Pk2b95MixYtSEtLY+vWrXzzzTckJSUVun9JFRZLoIMHD/r/wAc2vRYVb48ePXjzzTcZPnw4r7/+eolqQNu0acPatWvzbOvatStpaWns3buXmJgYZsyYwT333FPCq/M1zY4ePZq9e/fyr3/9K9/nBd3Xvn37cv311/Ob3/yGBg0asH//furXr1/oNfXp04fnnnuOe++9199EGfjvpqD7NGrUKK699lpSUlKIiYnJF1duzd35GjhwIFOmTOGmm27iyy+/pF69ejRp0oSrrrqK3//+9/73QD/88EP+/Oc/U79+ferWrcsXX3zBZZddxvTp0wu83wVdM8DWrVv5/PPP6d69O2+88QaXX375eV9DaYR4zV1cuZ5+/a5sNu72PdBa1SLp3VpNsiJSNaWkpJCZmUmvXr2IjIykWbNmhf5huu2220hMTCywxiXQtGnTGDt2LImJiaxYsYIJEyYAcN999/Hcc8/RoUOHPB0Cevfuzbfffltgh4q3336b+Ph4kpOTWbVqFSNGjKBBgwb07NmT+Ph4xo4dW6rrLiyWQL/97W/53e9+R4cOHfLUShYV79/+9jdeeeUVEhMTefXVV3nqqaeKjGPv3r0F1uI1adKEiRMn0rt3b5KSkujUqRODBg0q8fW1b9+ew4cP07RpU5o0yV8BUdB9bd++Pffffz9XXHEFSUlJ/s42hV3TU089xcKFC0lISKBTp07+IXNyJSYmEhkZSVJSEpMnTwZ8zfd169bl1ltvLfG1FOXpp58mNjaWzMxMEhMTGTVqFOAbLqZFixa0bNmS0aNH+3tX169fnwceeIAuXbrQpUsXJkyY4G9+fvbZZxk1ahQtW7bk0ksv5ZprrslXXmHX3Lp1a5555hnatm3LgQMH/E23Fc2q4vgtxel8UaRbelttSFsO9VuUWzlPfLiOpxdsBGBAYhOm3Nyx3MoKJ+np6aSmpgY7DCklPb+SWbNmDW3btg12GHlo+qrgmTt3Lps3byYtLa3U56hKz2/Hjh2kpqaydu1aIiJCo54pIyODAQMGsGrVqiL3K+i7b2bLnHOdyyqWkGyWPXLhxfBfz0HBMS43AAAWaElEQVS9ZsXvXErOOeYGNMkOUJOsiIiUUnEDJoeS6dOnc//99/PEE0+ETGJX2YRkcnc6ojq0vrpcy1j7/WE27zkCwAXVI0lVk6yIiEixRowYwYgRI4IdRpmLi4srttauoihlLqXAXrI/bduYmtUigxiNiIiIiI+Su1JQL1kRERGprJTclcK3Ow+xZa+vSfbC6pFc8ZNGQY5IRERExEfJXSkENsn2bacmWREREak8lNydo/xNshcFMRoRkfLTv39/srKyyMrK8o8PBr7hbkrbuzM9Pd0/sTz45jadPn36eccKcMsttxQ6EHGuqVOnsmPHmfnAR40alW9ctvKydOlS/1AnU6dO5e677y7xsWV5nwozZ84cJk6cCMDs2bPL/b5cffXVREdH5/u3tGXLFi677DJatmzJkCFDOH78eLnGEYqU3J2j1TsO8d2+owDUqRFFSquGQY5IRKR8zJ8/n+jo6HzJ3fk4O7m7/fbbK7Tn5NnJ3d///nfatWtXIWV37tyZp59+ulTHVsR9GjhwIOPHjwcqJrkbO3Ysr776ar7t48aNY8yYMWzcuJGYmBheeumlco0jFCm5O0dz1SQrIiHg8ccf9ycaY8aM4corrwR8E6LnzjwRFxfH3r17GT9+PJs2bSI5Odk/E0R2djaDBw+mTZs2DB06tMDZFZ5++mnatWtHYmIiN910ExkZGTz//PNMnjyZ5ORkFi9ezIMPPsikSZMASE1NZcyYMXTu3Jm2bdvy1VdfccMNN9CqVSv+8Ic/AL6BYuPj4/1lTJo0iQcffDBf2Q8//DBdunQhPj6e2267DeccM2fOZOnSpQwdOpTk5GR++OEHUlNTWbp0KQAzZswgISGB+Ph4xo0b5z9X7dq1uf/++0lKSqJbt27s2rWryHubk5PDrbfeSkJCgn8qNyi+xvP06dPExcXlmTu1VatW7Nq1K899Ovu+5j6PO+64g4SEBBITE3nnnXeKvKb333+fjh07kpSURJ8+fYAztYmfffYZc+bMYezYsSQnJ7Np0yY6djwzSP+GDRvyrJdWnz598g267JxjwYIFDB48GICRI0cye/bsfMdmZ2f773Hg9dauXZsxY8bQvn17+vTpw549e847zqooJMe5Ky++Jtkz/8enXrIiUiYezD9/a9md+2CBm1NSUvjrX/9KWloaS5cu5dixY5w4cYLFixfTq1evPPtOnDiRVatW+ef6TE9PZ/ny5axevZqLLrqInj178u9//zvfdGUTJ05ky5Yt1KhRg6ysLKKjo7n99tupXbs29913HwCffPJJnmOqV6/O0qVLeeqppxg0aBDLli2jfv36XHrppYwZM6bEl3333Xf7pzsbPnw4c+fOZfDgwUyZMoVJkybRuXPeyQB27NjBuHHjWLZsGTExMfTr14/Zs2dz3XXXceTIEbp168af/vQnfvvb3/Liiy/6k82CPPPMM5gZK1euZO3atfTr14/169cXG3NERASDBg1i1qxZ3HrrrXz55ZdccsklNG7cOM9+Z99XgEceeYS6deuycuVKAA4cOFDoNfXs2ZPRo0ezaNEimjdvzv79+/Ocv0ePHgwcOJABAwb4k6x69eqxYsUKkpOTeeWVVwqcNuzxxx/n9ddfz7e9V69eJa6x3LdvH9HR0URF+dKT2NhYtm/fnm+/Rx55hHr16uW5XoAjR47QuXNnJk+ezMMPP8xDDz3ElClTSlR2KFHN3TlYtf0Q2/b/AECdmlFcriZZEamiOnXqxLJlyzh06BA1atSge/fufP311yxevLhEk9x37dqV2NhYIiIiSE5OJiMjI98+ufPPvvbaa/4/1sUZOHAgAAkJCbRv354mTZpQo0YNWrRowbZt20p8fQsXLuSyyy4jISGBBQsWsHr16iL3/+qrr0hNTaVRo0ZERUUxdOhQFi1aBPgSztwat06dOhV4rYE+/fRThg0bBkCbNm245JJLSpTcAQwZMsQ/T+2bb77JkCFD8u1T0H39+OOPGT16tH+fmJiYQq/piy++oFevXjRv3hzAP6dqUUaNGsUrr7zCqVOneOutt7j55pvz7TN27FhWrFiR76e0TdFF+fjjj7nrrrv86zExMYAvQc69Z8OGDePTTz8t87KrAiV352BuQK1dv3Y/pkaUmmRFpGqqVq0azZs3Z+rUqfTo0YOUlBQWL17Mxo0bSzTnbY0aNfzLkZGRnDx5Mt8+8+bN46677uLrr7+mS5cuBe5T2HkjIiLylBEREcHJkyeJiori9OnT/u05OTn5zpGTk8Odd97JzJkzWblyJaNHjy5wv5KqVq0aZgYUfq1lpXv37mzcuJE9e/Ywe/Zsbrjhhnz7lOa+nq+f//znvPfee8ydO5dOnTrRoEGDfPs8/vjjJCcn5/s5l/lyGzRoQFZWlv+aMjMzadq0aanjzn1u4UbNsufgvZXf+5c1l6yIlJlCmk7LW0pKCpMmTeLll18mISHB/77b2X8Q69Spw+HDh8/p3KdPn2bbtm307t2byy+/nDfffJPs7Gzq1KnDoUOHSh1z48aN2b17N/v27aN27drMnTuXq6/OO91kbiLXsGFDsrOzmTlzpr95sbBr6dq1K2lpaezdu5eYmBhmzJjBPffcU2Qss2bNYsmSJfz5z3/Osz0lJYXXX3+dK6+8kvXr17N161Zat27N559/Xuz1mRnXX389v/nNb2jbtm2+JKqw+9q3b19efPFFf8eXAwcOFHpN3bp1484772TLli3+Ztmza+/Ovk81a9bkqquu4o477ii0g8PYsWP972SWlpnRu3dvZs6cyU033cS0adMYNGhQvv369u3LM888w5NPPum/3piYGE6fPu0/9o033sj3qkC4UM1dCR07eYqt+329ZCMjjJ4t1SQrIlVbSkoKO3fupHv37jRu3JgaNWoU2CTboEEDevbsSXx8fIn/eJ86dYphw4b5OxWkpaURHR3Ntddey6xZs/wdKs5VtWrVmDBhAl27dqVv3760adMm3z7R0dGMHj2a+Ph4rrrqKrp06eL/7JZbbuH222/3d6jI1aRJEyZOnEjv3r1JSkqiU6dOBSYVgTZt2kTdunXzbb/zzjs5ffo0CQkJDBkyhKlTp+aphSzOkCFDeO211wpski3svv7hD38gKyuL+Ph4kpKSWLhwYaHX1KhRI1544QVuuOEGkpKSCiznpptu4vHHH6dDhw5s2rQJgKFDhxIREUG/fv1KfC1FSUlJ4cYbb+STTz4hNjaWDz74AIDHHnuMJ554gpYtW7Jv3z5+9atf5Tv2D3/4AwcOHMhzvQAXXnghS5YsIT4+ngULFvjfuww3VlAPp6qudevWbt26dWV6zv1HjtPxkY8AiL6gGismlM0/bskvPT2d1NTUYIchpaTnVzJr1qwpUfNnRTp8+HC+3otSuGHDhjF58mQaNaocsxSV9/ObNGkSBw8e5JFHHim3Ms5X7dq1yc7ODnYYRSrou29my5xznQs55JypWbaEsnPOvNNwYXXdNhGRcPfaa68FO4QKc/3117Np0yYWLFgQ7FCkBJSllFD2sTPJXe0aum0iIhI+Zs2aFewQSqSy19pVFL1zV0JHjgckdzWV3ImIiEjlpOSuhPI0y6rmTkTKQCi+8ywihauo77ySuxIKbJato+RORM5TzZo12bdvnxI8kTDhnGPfvn3UrFmz3MtSllJCgcndhTU0eLGInJ/Y2FgyMzMr1dyXOTk5FfKHR8qHnl/lV7NmTWJjY8u9nCqT3JnZ1cBTQCTwd+fcxIos/8gxNcuKSNnJnSGiMklPT6dDhw7BDkNKSc9PclWJZlkziwSeAa4B2gH/ZWbtKjIGNcuKiIhIVVAlkjugK7DRObfZOXcceBMoeujwMqYOFSIiIlIVVJXkrimwLWA909tWYQKHQlFyJyIiIpVVyGQpZnYbcJu3eszMVpVXWcMfg+HldXIBaAjsDXYQUmp6flWXnl3VpudXdbUuy5NVleRuO9AsYD3W2+bnnHsBeAHAzJaW5RxtUrH0/Ko2Pb+qS8+uatPzq7rMbGlZnq+qNMt+BbQys+ZmVh24CZgT5JhEREREKp0qUXPnnDtpZncDH+AbCuVl59zqIIclIiIiUulUieQOwDk3H5hfwt1fKM9YpNzp+VVten5Vl55d1abnV3WV6bMzTX0jIiIiEjqqyjt3IiIiIlICIZfcmdnVZrbOzDaa2fhgxyN5mVkzM1toZt+a2Woz+7W3vb6ZfWRmG7zfMd52M7Onvef5jZl1DO4VCPhmjTGz5WY211tvbmZfes/pLa/jE2ZWw1vf6H0eF8y4Bcws2sxmmtlaM1tjZt31/asazGyM99/NVWY2w8xq6rtXeZnZy2a2O3BottJ818xspLf/BjMbWZKyQyq5qwzTlEmxTgL/zznXDugG3OU9o/HAJ865VsAn3jr4nmUr7+c24LmKD1kK8GtgTcD6Y8Bk51xL4ADwK2/7r4AD3vbJ3n4SXE8B7zvn2gBJ+J6jvn+VnJk1BdKAzs65eHydC29C373KbCpw9Vnbzum7Zmb1gT8Cl+GbreuPuQlhUUIquaMSTFMmRXPO7XTOfe0tH8b3h6Upvuc0zdttGnCdtzwImO58vgCizaxJBYctAcwsFvgZ8Hdv3YArgZneLmc/v9znOhPo4+0vQWBm9YBewEsAzrnjzrks9P2rKqKAWmYWBVwA7ETfvUrLObcI2H/W5nP9rl0FfOSc2++cOwB8RP6EMZ9QS+6CPk2ZlJzXTNAB+BJo7Jzb6X30PdDYW9YzrXyeBH4LnPbWGwBZzrncOfoCn5H/+XmfH/T2l+BoDuwBXvGa1f9uZhei71+l55zbDkwCtuJL6g4Cy9B3r6o51+9aqb6DoZbcSRVhZrWBd4B7nXOHAj9zvi7c6sZdCZnZAGC3c25ZsGORUokCOgLPOec6AEc40ywE6PtXWXlNcYPwJegXARdSghocqbzK87sWasldsdOUSfCZWTV8id3rzrn/8zbvym3u8X7v9rbrmVYuPYGBZpaB77WHK/G9wxXtNRVB3mfkf37e5/WAfRUZsOSRCWQ657701mfiS/b0/av8fgpscc7tcc6dAP4P3/dR372q5Vy/a6X6DoZacqdpyio5752Pl4A1zrknAj6aA+T2AhoJvBuwfYTXk6gbcDCgSlsqmHPud865WOdcHL7v1wLn3FBgITDY2+3s55f7XAd7+6tWKEicc98D28wsd5LyPsC36PtXFWwFupnZBd5/R3Ofnb57Vcu5ftc+APqZWYxXe9vP21akkBvE2Mz643snKHeasj8FOSQJYGaXA4uBlZx5Z+v3+N67exu4GPgO+IVzbr/3H7Ep+JofjgK3OufKdIJlKR0zSwXuc84NMLMW+Gry6gPLgWHOuWNmVhN4Fd+7lfuBm5xzm4MVs4CZJePrDFMd2Azciu9/9PX9q+TM7CFgCL5RB5YDo/C9f6XvXiVkZjOAVKAhsAtfr9fZnON3zcx+ie/vJMCfnHOvFFt2qCV3IiIiIuEs1JplRURERMKakjsRERGREKLkTkRERCSEKLkTERERCSFK7kRERERCiJI7ERERkRCi5E4kDJnZKTNbYWarzOyfZhZ9jsenm1lnb3n+uR5fyDlnmNk3ZjbmfM9VRBnRZnZnwHqqmc0tr/KKiMNfrpndYmZTvOWpZja46KPLJZ5aZvYvM4v01h83s9Vm9ngpzpXsjTeauz7AzB4uy3hFpGhK7kTC0w/OuWTnXDy+AU7vKu2JnHP9nXNZ5xOMmf0Y6OKcS3TOTT7rs6hCDiuNaODOYveqgkpznwKO+SXwf865U976bUCic25sKUJJBvoHrM8DrjWzC0pxLhEpBSV3IvI5vlHuMbOuZva5mS03s89yp6nyanbeNLM1ZjYLqJV7sJllmFlDM4szs1UB2+8zswe95TQz+9armXuzgBg+BJp6tYkpXs3gk2a2FPi1mfXxYlppZi+bWY2Asv/sHbfUzDqa2QdmtsnMbi+gnInApd7+ubVStc1sppmtNbPXvZHiMbNOXm3WMu+cTQJPZGaRZrbFmy4o2qsN7eV9tsjMWpnZhV68S7z4B5XgefzUu5b1ZjbAO19NM3vFu/7lZtbb236Lmc0xswXAJyUpz6s1XGxmc/BNXwUwFG8aJG97bWCZmQ0xs0Zm9o6ZfeX99PT2y1eW+aZ9fBgY4t3jId6UV+nAgBJcu4iUgbL8P2IRqWK8Zrg++Ob7BVgLpDjnTprZT4FHgZ8DdwBHnXNtzSwR+PocixoPNPemRSqoCXcgMNc5l+zFBVDdOdfZfNMobQD6OOfWm9l0L54nvWO3OueSzWwyMBXfZOo1gVXA8wXEER9QTiq+6ZnaAzuAfwM9zexL4G/AIOfcHjMbAvwJXw0XAM65U2a2DmgHNPfuSYp3bDPn3AYzexTfnJ6/9K57iZl9XMy9igO6ApcCC82sJb6aVeecSzCzNsCHZvYTb/+O+GrZ9hdWnnPuyFlldPTuwxYvIWvhnMvwrmugmWUH3KM3gMnOuU/N7GJ881q2Be4/uyzgY2AC0Nk5d3dAeUuBFHzTLolIOVNyJxKeapnZCnw1dmuAj7zt9YBpZtYKcEA1b3sv4GkA59w3ZvbNOZb3DfC6mc3GN7diSbzl/W4NbHHOrffWp+FLdnKTuzne75VAbefcYeCwmR0zs+gSNBkvcc5lAnj3JA7IAuKBj7xEMxLYWcCxi/Hdm+bAn4HRwL+Ar7zP+wEDzew+b70mvjkli/K2c+40sMHMNgNtgMvxJZs459aa2XdAbnL3kXNufzHlrSngmrd4yw296y3MT4F23n0AqGtmtc/x2nYDFxVRhoiUISV3IuHpB6+26wJ8NTF34UveHgEWOueuN7M4fM1pJXWSvK961AxY/hm+JOha4H4zS3DOnSzmfGfXNhXmmPf7dMBy7npJ/hsXeMwp7xgDVjvnuhdz7CJ8tYgX4auxGotvovDF3ucG/Nw5ty7wIDNrXMQ5z57wu7gJwAPvU4HlFXPMD+R9VmeLALo553ICN3rN1wVd22UFnKOmV46IVAC9cycSxpxzR4E04P+Z7+X6esB27+NbAnZdBNwMYGbxQGIBp9sF/MjMGnjvxOW+LxaBr5lyITDOK6P2OYS5DojzmicBhuOrHSuNw0CdEpbZyMy6A5hZNTNrX8B+S4AewGkv+VkB/De++wW+xPmegPf4OpSg7BvNLMLMLgVaeLEsxvdeHF5z7MXe9rMVWJ6ZNTWzTwoqzDl3AIj0mr8L8iFwT+6KmSUXc20F3eOf4GsmF5EKoOROJMw555bjazb9L+AvwJ/NbDl5a72ew9fxYA2+F+aXFXCeE95nS/A18671PooEXjOzlcBy4Olz6V3rJU23Av/wznGa/O/SlfRc+4B/m28ImEKH+XDOHQcGA4+Z2X/wJW09CtjvGLAN+MLbtBhfYrPSW38EX9P2N2a22lsvzlZ89/A94Hbv+p8FIrzrfwu4xSv7bIWV1wRfzWphPsTX9FuQNKCz+TrDfAvkdlQprKyF+JpxV3jvKgL0xtdrVkQqgPk6MomISKgys7vxdTyZU8jnHYExzrnh5VB2Y+AN51yfsj63iBRMyZ2IiGBmvwSmBYx1V1bn7QKccM6tKMvzikjhlNyJiIiIhBC9cyciIiISQpTciYiIiIQQJXciIiIiIUTJnYiIiEgIUXInIiIiEkL+P5uZLtGxuDQRAAAAAElFTkSuQmCC\n",
            "text/plain": [
              "<Figure size 720x432 with 1 Axes>"
            ]
          },
          "metadata": {
            "tags": [],
            "needs_background": "light"
          }
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "vY-LcPAG03w4",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 34
        },
        "outputId": "fbe88e69-a49f-44f5-b439-34ff6f3ec91b"
      },
      "source": [
        "#Values used : \n",
        "\n",
        "[u,u_stim,q,pe,B,re]"
      ],
      "execution_count": 67,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "[500, 200, 200, 3000, 1.125, 1000]"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 67
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "eKKHccTUwgzX",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        ""
      ],
      "execution_count": null,
      "outputs": []
    }
  ]
}